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Abstract 

In this paper, we show that there exists an essential difference of boundary effects between 
Bose and Fermi systems both for Dirichlet and Neumann boundary conditions: at low tem- 
peratures and high densities the influence of the boundary on the Bose system depends on 
the temperature but is independent of the density, but for the Fermi case the influence of the 
boundary is independent of the temperature but depends on the density, after omitting the 
negligible high-order corrections. We also show that at high temperatures and low densities 
the difference of the influence of the boundary between Bose and Fermi systems appears in 
the next-to-leading order boundary contribution, and the leading boundary contribution is in- 
dependent of the density. Moreover, for calculating the boundary effects at high temperatures 
and low densities, since the existence of the boundary modification causes the standard virial 
expansion to be invalid, we introduce a modified virial expansion. 

PACS numbers: 05.30.-d, 05.70.Ce, 65.80.-|-n, 82.60.Qr 

1 Introduction 

In small-size and low-temperature systems, the mean thermal wavelength of particles is comparable 
with the size of the system, and the boundary effect becomes one of the most important quantum 
effects. The influence of boundaries is widely discussed. Based on Maxwell-Boltzmann statistics, 
the boundary effect on the classical ideal gas confined in a narrow box or in spherical and cylindrical 
geometries is studied theoretically in 1 , 2 . For quantum gases, the influence of the boundary on 
Bose gases is discussed in 0; the shape and topology dependence of the boundary effect on ideal 
quantum gases confined in irregular containers is discussed in jl]; some exact and approximative 
solutions for quantum gases in finite-size containers are provided in a general result for the 
boundary modification on ideal quantum gases in arbitrary dimensions is provided in |6|. In 
Bose systems, the choice of the boundary condition may play an important role. Especially, for 
Bose systems it may display a phase transition — the Bose-Einstein condensation. Systematic 
discussions on this subject are given in 013121 EH- Moreover, some experimental studies (e.g., 
[ni 112 Ennuis') show that the influence of the boundary may be remarkable in low-temperature 
quantum gases confined in small volumes. 

In the present paper we pay special attention to two-dimensional systems. Many novel proper- 
ties of the gases absorbed within a bundle of carbon nanotubes are reported and, especially, 
under appropriate thermodynamic conditions gases adsorbed within the nanotubes and on the 
external surface of the bundle will display two-dimensional behaviour ^|E1EI- Moreover, we 
also provide some three-dimensional results. 
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Different boundary conditions lead to different boundary effects. Two kinds of boundary con- 
ditions are considered: the Dirichlet boundary condition and the Neumann boundary condition. 

The primary motivation of this work is to compare the boundary effects on Bose and Fermi 
systems. The analysis given by [5] shows that the ideal quantum gas can be used as an efficient tool 
for studying the boundary effects in many quantum systems, including electron gases, interacting 
Bose gases, etc, since in an interacting system there are two kinds of effects: one arises from the 
classical interaction, e.g., the van der Waals interaction, and the other is the quantum effect. In 
ideal quantum gases though there are no classical interactions the quantum effects, such as the 
exchange interaction, remain. The boundary effect is, however, a kind of quantum effect, which 
becomes important when the mean thermal wavelength of particles is comparable with the size of 
the system, so the boundary effects on ideal and interacting gases are similar and the conclusions 
of the ideal gases are of widespread applicability. 

In finite volumes (areas), energy levels are discrete, and the spectra of the particles are shape 
dependent and sensitive to the topology To take boundary effects into account, strictly speak- 
ing, we need to perform the summation over all possible states directly, but that is in general 
hardly to be solved. In 0] and [S], some methods are developed to perform the summation. Based 
on a mathematical work given by M. Kac [20], a method for calculating the boundary modification 
to the ideal quantum gases in irregular-shaped containers is provided in ^ ; some thermodynamic 
quantities for two- and three-dimensional systems are calculated in [3] and some of the results are 
exact solutions. Different boundary conditions lead to different boundary effects. The equation of 
state for ideal quantum gases with boundary modifications is provided in 0] and |HI . The boundary 
condition used in these two papers is the Dirichlet boundary condition. In the following, we will 
consider the boundary effects corresponding to two kinds of boundary conditions: the Dirichlet 
boundary condition and the Neumann boundary condition. It is straightforward to obtain the 
equation of state with the Neumann boundary condition by use of the method given in 0] and 
0. The equation of state of Fermi and Bose gases in confined space with Dirichlet and Neumann 
boundary conditions can be expressed as 

— = Y^A,K{z), (1) 

CT 

N = Y,AaK-iiz), (2) 

(T 

where h^r (z) equals the Bose-Einstein integral (z) or the Fermi-Dirac integral /o- (z) in the Bose 
or the Fermi case, respectively: 
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In these equations and the following, the upper sign stands for the Dirichlet boundary condition 
and the lower sign for the Neumann boundary condition. In two dimensions, — S is the area, 
CT = 2, 3/2, 1, and A2 = gS/X^, A3/2 = T (1/4) gL/X, Ai = gx/6, where A = h/V2TrmkT is the 
mean thermal wavelength, L is the perimeter of the boundary, x is the Euler-Poincarc character- 
istic number (the Euler-Poincare characteristic number reflects the connectivity — a topological 
property — of a two-dimensional system), and g denotes the number of internal degrees of free- 
dom, the degeneracy number of spin states (for bosons we take g = 1). For the case of the gas 
confined in a three-dimensional box, = ^ is the volume, ct = 5/2,2,3/2,1, and = gV/X^, 
^2 = T (1/4) gS/X^, A3/2 = (1/16) gi/A, Ai = ^ (1/8) g, where S is the area of the surface of the 
box, L — 4:{Lx + Ly + Lz) is the total length of the sides of the box. 
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The difference between the Dirichlet case and Neumann case is that the leading contributions 
of these two kinds of boundaries have opposite signs: the contribution of the Dirichlet boundary 
is negative, and the contribution of the Neumann boundary is positive. 

The standard virial expansion approach cannot be used due to the existence of the boundary 
modification. For analysing the boundary effect in high-temperature and low-density systems, we 
introduce a modified virial expansion. 

In this paper, we compare the boundary effects on the Bose system and on the Fermi system. 
At low temperatures and high densities, our results will show that the boundary effects on Bose 
and Fermi systems are different essentially: after omitting the negligible high-order corrections, 
for Bose systems the influence of the boundary depends on the temperature but is independent 
of the density; however, contrary to the Bose case, the influence of the boundary on Fermi sys- 
tems depends on the density but is independent of the temperature. At high temperatures and 
low densities, the difference of the boundary effects between Bose and Fermi systems appears in 
the next-to-leading boundary contributions: the next-to-leading boundary contributions to Bose 
systems and to Fermi systems have opposite signs; they are functions of both temperature and 
density, though the leading boundary contribution is independent of the density. 

In Section 2, we compare and discuss the difference of boundary effects between Bose and 
Fermi systems both for Dirichlet and Neumann boundary conditions. In Section 3, by introducing 
a modified virial expansion, we calculate the boundary effects on Bose and Fermi systems at high 
temperatures and low densities. The conclusions are summarized in Section 4. 

2 The difference between Bose and Fermi systems at low 
temperatures and high densities 

In this section, we compare the boundary effects on the two-dimensional Bose and Fermi systems 
at low temperatures and high densities. We show that the boundary effects on Bose and Fermi 
systems are different essentially, and discuss the reason of such a difference. 

We will take the specific heat, a quantity directly accessible to experimental measurement, as 
an example, and our main conclusion holds also for other thermodynamic quantities. 

The Bose case. For Bose cases, from the equation of state, Eqs. Q and (jSjl, we can ob- 
tain the specific heat at low temperatures and high densities. From Eq. (O one can determine 
a = 1 ^ {^/n/A) (LX/S) e^^^ 1"^ e~"'^^ approximately, where a = —[ij (kT), with fi the chemical 

potential, and n — N/S is the mean number density, where N is the expectation value of the 
particle number in the grand canonical ensemble. Substituting a into Eq. (^), we obtain the grand 
potential. Starting from the grand potential we can calculate the specific heat. 
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(4) 

In this section, we always drop out of the topological terms which are often negligible. The terms 
proportional to S/X'^ in Eq. Q are the specific heat in free space, and the remaining terms are 
the boundary modifications. 

The percentage of the boundary modification to specific heat is 



Cv ~ Cy „ , ,„ i A 
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1 - 0.226 (nA^) e 
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We can see that the boundary effect on a Bose system, after omitting the negligible high-order 
corrections, is independent of the density; it is proportional to X/\/S, which is the ratio between 
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the thermal wavelength and the linear size of the container, and L/^/S^ which reflects to some 
extent the information of the shape of the boundary. We can see that as the temperature falls, the 
contribution from the boundary will become more and more important. 

The above result is obtained in low temperatures. It should be emphasized that this result 
will become invalid when the fugacity z ^ 1, since the result which we based on given by Refs. 
mi and |S] is invalid when z 1. If the system can display Bose-Einstein condensation, the case 
z ^ 1 corresponds to a non-zero temperature T^. That is to say, if the Bose-Einstein condensation 
occurs, our result will be invalid at a certain non-zero temperature. However, in the above case 
(the two-dimensional ideal Bose gas with Dirichlet and Neumann boundaries) there is no Bose- 
Einstein condensation [7]. This means that only when T — > the fugacity z ^ 1. Nevertheless, 
the modification to the equation of state, Eqs. (Q) and will be valid only when the boundary 
contribution is small [201 ■ Consequently, for the Bose case there exists a lower limit on temperatures 
in the range of the applicability of the above result. This lowest temperature can be roughly 
estimated: the boundary modification should be less than the leading contribution. From Eq. 
we can obtain this lower limit for a square two-dimensional system directly: the above result is 
valid only when the temperature T > 0.057/i^/ (mka^), where a is the side length of the system. 
For illustration, such a temperature for ^■^Na with size 50nm is about 2 x 10~^K. Moreover, the 
boundary influence on the condensation is a very important and interesting subject, which has 
been considered in [7ll5ll^[Tn|. 

The Fermi case. For Fermi cases, the specific heat is 



Cv_ 
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The ratio of the boundary modification to the specific heat in free space is 



Cv 
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where 



2\J ny/Sy/S 



Unlike the Bose case, after omitting the contribution suppressed by (kT/e^Y, the ratio {Cv ~ 
(^v)/^v TL/VS{\/^/S)yJ S/N , is independent of the temperature but is determined by such 
three factors: the number density N/S (or, the mean space between particles yJS/N), the shape, 
which is described by L/^/S to some extent, and the linear size v^. In other words, the boundary 
effects on two-dimensional Fermi systems depend on the density, the size, and the geometrical 
properties of the system, but is independent of the temperature. 

The difference between Bose and Fermi systems. In a word, in addition to the size and shape, 
the boundary effect on Bose systems depends almost only on the temperature, but that on Fermi 
systems depends almost only on the density. Such a difference is quite essential. In the following 
we will analyse the reason. 

Essentially, the influence of the boundary is determined by the thermal wavelength of the 
particles and the size of the system. The lower the energy of a particle, the longer is the wavelength 
and then the stronger is the boundary effect. Therefore, the particles in lower energy levels will 
be influenced more strongly than those in higher energy levels. 

First, we discuss the influence of the temperature on the boundary effects. 

For a Bose system, when the temperature falls, the particles will tend to occupy lower energy 
levels, so the boundary effect will become stronger. As a result, the boundary effect will be a 
function of the temperature. 

For a Fermi system, there exists the Fermi energy. At low temperatures almost all the particles 
are in the states below the Fermi energy. The magnitude of the influence of the boundary on the 
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particles is determined by the wavelength: the shorter the wavelength the weaker the boundary 
influence is. Generally speaking, the energy of the particles near the Fermi surface is relatively high, 
so the wavelength is short. Therefore in a low-temperature system the contribution of the boundary 
influence mainly comes from the particles far below the Fermi surface. At very low temperatures, a 
small change in the temperature only influences the particles near the Fermi surface; however, due 
to their short wavelengths, the boundary contribution from such particles is very small compared 
with the whole boundary effect. Therefore, in low-temperature Fermi systems the boundary effect 
is almost independent of the temperature. 

Next, we discuss the influence of the density on the boundary effects. 

For a Bose system with fixed size and shape, the change in the density will change the total 
number of particles, but almost has no influence on the relative distribution of particles. What 
we consider is the ratio between the boundary contribution and the result of free space, i.e., the 
percentage of the boundary contribution, which mainly depends on the distribution of the particles. 
Therefore, such a ratio will not change after changes in density. 

For a Fermi system with fixed size and shape, the ratio of the boundary contribution to the 
free space result will get smaller as the density increases. When the volume is given, as the 
density increases, the number of particles increases. Roughly speaking, the energy of the newly 
added particle will be higher than the original Fermi energy, i.e., the wavelength of the newly added 
particle will be shorter than the original mean wavelength. As mentioned above, the boundary effect 
is mainly determined by the magnitude of the wavelength of the low-energy particles. Therefore 
with the density increasing the proportion of the influence of the boundary with respect to the 
result in free space will be suppressed. 

Particularly, it should be emphasized that in the Fermi case the influences on different thermo- 
dynamic quantities are different: for example, for the Dirichlet boundary condition, the influence 
on the specific heat is negative, but on the Fermi energy is positive. For illustrating this, we give 
the chemical potential and Fermi energy in a two-dimensional system, 
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(8) 



eF = e%il±S), (9) 

where — (?i^/2to) {inN/ gS) is the Fermi energy in two-dimensional free space, and [ep — e^) /e^ 
zLS. Note that this boundary modification to the case of the Dirichlet boundary condition is posi- 
tive and to the case of the Neumann boundary condition is negative; however, for the specific heat, 
Eq. O, the modification to the Dirichlet case is negative, but to the Neumann case is positive. 

Moreover, the boundary effect on the chemical potential deserves special attention. For illus- 
trating this, we give the three-dimensional chemical potential and Fermi energy. 



/"3 — £_F3 



(1± A) {-^ 



F3 



(10) 



and 

£f 3 = 43 (1 ± A) , (11) 

where e^^^g — {h^ /2rn) [Q-k'^N / gV^^^ is the Fermi energy in three-dimensional free space and 

^_l/7r5N^i/3 S 1 fV^^'"^ 



4 V 6 / \n 

The second-order term of the two-dimensional chemical potential (|SJ) only includes the boundary 
contribution, while the second-order term of the three-dimensional chemical potential (|10|l is inde- 
pendent of the boundary, i.e., there is no second-order boundary correction to the three-dimensional 
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chemical potential. This leads to such a result: in two dimensions, the second-order contribution 
to the ratio (/^ — /^") = ±d 1 — (7r^/24) (kT /e'pY + ' ' ■ is negative, but in three dimensions 

the second-order contribution to the ratio (/ua — ^3) //i-j = ±A 1 -t- (7r^/12) {kT/ep^Y + ' ' ' is 
positive. The reason is that in the result of the standard statistical mechanics, which is under the 
thermodynamic limit approximation, the two-dimensional chemical potential has no second-order 
term [2], so the second-order contribution in Eq. © only comes from the boundary effect. In 
three-dimensional free space, however, the chemical potential does have the second-order term, 
but the three-dimensional boundary correction, which behaves like the two-dimensional chemical 
potential, has no second-order term. 

3 Boundary effects at high temperatures and low densities: 
a modified virial expansion 

The existence of the boundary modification causes the standard virial expansion to be invalid. For 
analysing the boundary effect at high temperatures and low densities, we introduce a modified 
virial expansion approach. 

In statistical mechanics, in the limit T — > 00, hence A — > 0, all the thermodynamic quantities can 
be expanded with respect to nX^ since nX^ — NX^/S <C 1, which is the so-called virial expansion. 

The main idea of the virial expansion in two dimensions is to suppose that the equation of state can 

00 

be expanded as PS/{NkT) — a; (nX^) , where the coefficients a/ are referred to as the virial 

1=1 

coefficients. Nevertheless, since Eqs. ^ and (|2l contain the terms which describe the influence of 
the boundary shape and topology, the expansion parameter of the virial expansion nX^ (or NX^ /S) 
is not the unique expansion parameter, so the thermodynamic quantities cannot be expressed as 
the series of nX^. In other words, it is impossible to perform the standard virial expansion. 

For calculating the thermodynamic quantities at high temperatures and low densities, we intro- 
duce a new approach which is, in fact, a modified virial expansion. We suppose that the equation 
of state can be expressed as a series, 

PS , ,2n('-1)/2 



NkT 



^a,+l)/,(nA^)^'-^^/^ (12) 



where a(i_|_i)/2 are the modified virial coefficients. 
Using Eqs. ((U and (PJ, we obtain 

nA'= z + 77— + ••• Tt^ ^ + r;— + ■•• + r;z2 + ...). (13) 



2 J ' A S \ 'V2 J & S 

In this equation and following, for the Bose case 77 = +1 and for the Fermi case 77 = —1. In this 
equation, the expansion parameters include NX^/S, LX/S, A"^/ (LS) and (A^/S')x/6. Therefore 
we write z in the form 

- = -^'EliW2(-A^)'"'^', (14) 

where fo(;+i)/2 are the coefficients of the expansion. Then we can obtain an equation of A. By 
equating the coefficients of each power of A, we can obtain the coefficients &(;+i)/2 and then obtain 
the modified virial coefficients a(;+i)/2: 

ai = 1, a3/2 = 0, 02 = -0.25r,, a^/^ = T0.037r/^ A , ^3 ^ 0.028 - O.OOSr;!^, • ' ' • (15) 

V JV V "J iv D 

One can see that a(;+i)/2 a-re related to the geometry of the system, and if the thermodynamic 
limit approximation is taken, they will return to the common virial coefficients. By a(;^i)/2 7 we can 
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express the expansions for the thermodynamic quantities. As an example, we give the expansions 
of chemical potential and specific heat: 



= In (nA2)± 0.25^ + 0.031(^1 -0.167^ 



.LA 



kT 



LX 



S 



1± 0.146^ + 0.011 



S 



LAV 



0.042 {nX^y 



LX 

1 ± 0.100— + ■ 



+ ■ 



Nk 



LX 

1 ± 0.063— + 



0.028 (nA^ 



1 + 0.188 



0.25r?nA^ 
LX 



LA /LA 
+0.110— + 0.021 i — 



In comparison with the virial expansion, one can see that the boundary effects are reflected 
mainly by the terms which are proportional to LX/ S and its powers. 

The result shows that at high temperatures and low densities, both for Bose and Fermi cases, 
the leading contribution of the boundary to the chemical potential and the specific heat is a 
function of temperature and geometrical property of the system, but is independent of the density 
of the system. In the next-to-leading order boundary contribution, the difference between Bose 
and Fermi systems appears: the modifications to Bose systems and to Fermi systems have opposite 
signs. Moreover, the boundary modifications to different thermodynamic quantities are different. 
One can also verify that the properties we have drawn in this section are also satisfied by the other 
thermodynamic quantities. 

Note that the same procedures can be used to the three-dimensional case directly by expressing 
the equation of state and the fugacity z as: 



PV 
NkT 



1=1 



{nX^f-'^" 



and 



nA 



oo 
1=1 



{nX^f-'^'\ 



(16) 



4 Conclusion 

In conclusion, we show the essential difference of boundary effects on Bose and Fermi systems in 
two dimensions: at low temperatures and high densities, for Bose systems, boundary effects are 
almost independent of the density but depend on the temperature; while for Fermi systems, the 
effects are almost independent of the temperature but depend on the density. At high temperatures 
and low densities, the difference of boundary effects between Bose and Fermi systems appears in 
the next-to-leading order boundary contributions: the next-to-leading order contributions to Bose 
systems and to Fermi systems have opposite signs; moreover, the leading boundary contribution is 
independent of the density. 

Strictly speaking, for a finite-size system one should use the canonical formalism, but it is 
shown in |S] that the result based on the grand canonical formalism for finite-size systems is valid 
for most realistic cases. A detailed discussion for the canonical ensemble formalism is given by 
|22[I23| . In this paper attention is concentrated on the difference of boundary effects between Bose 
and Fermi systems, and we have not considered the Bose-Einstein condensation since there is no 
corresponding phase transition in the Fermi system at low temperatures. Of course, the boundary 
effect on the Bose-Einstein condensation is very important. Systematic discussions on this subject 
can be found in [7| HI |51 ^]. Note that in this paper though our attention focuses mainly on 
the two-dimensional cases, it is easy to check that most of our conclusion is also valid for three 
dimensions. 
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